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DESIGN OF TWO-DIMENSIONAL CHANNELS WITH PRESCRIBED VELOCITY DISTRIBUTIONS 

ALONG THE CHANNEL WALLS 1 

By John D. Stanitz 


SUMMARY 

A general method of design is developed for two-dimensional 
unbranched channels with prescribed velocities as a function of 
arc length along the channel walls. The method is developed 
for both compressible and incompressible, irrotational, non- 
viscous flow and applies to the design of elbows, difusers, noz- 
zles, and so forth. Two types of compressible flow are con- 
sidered: the general type, with the ratio of specific heats y equal 
to 1.4, for example, and the linearized type, in which y is —1.0. 
Two methods of solution are used: In part I solutions are 
obtained by relaxation methods; in part II solutions are 
obtained by a Green’ s function. 

Five numerical examples are given in part I including three 
elbow designs with the same prescribed velocity ass a function of 
arc length along the channel walls bid with incompressible, 
linearized compressible, and compressible flow. It is con- 
cluded that if a nonviscous gas with arbitrary y (1 . 4 , for ex- 
ample) were to flow through a channel designed for linearized 
compressible flow (7= —1.0), the resulting velocity distribution 
along the channel walls would be nearly the velocity distribution 
prescribed for the linearized compressible flow. 

One numerical example is presented in part II for an acceler- 
ating elbow with linearized compressible flow, and the time 
required for the solution by a Green’s function in part II was 
considerably less than the time required for the same solution 
by relaxation methods in part I. 

INTRODUCTION 

There are two general types of theoretical problem in 
two-dimensional fluid motion: (1) the direct problem, in 
which the distribution of velocity is determined for a pre- 
scribed shape of boundary, and (2) the inverse problem, in 
which the shape of boundary is determined for a prescribed 
distribution of velocity along the boundary. The direct 
problem is an anlaysis problem; the inverse problem is a 
design problem. This report is concerned with the inverse, 
or design, problem for two-dimensional, irrotational flow in 
unbranched channels with prescribed velocities as a function 
of arc length along the channel walls. 

The design of channels with prescribed velocities is impor- 
tant because: (1) Boundary-layer separation losses can be 
avoided by prescribed velocities that do not decelerate rap- 
idly enough to cause separation, (2) shock losses in com- | 


pressible flow and cavitation in incompressible flow can be 
avoided by prescribed velocities that do not exceed certain 
maximum values dictated by these phenomena, and (3) for 
compressible flow the desired flow rate can be assured- by 
prescribed velocities that do not result in “choke flow” 
conditions. 

Several methods of channel design have been developed 
for particular application (refs. 1 and 2, for example). In 
reference 1 a design method is developed for accelerating 
elbows in which the velocity increases monotonicafly along 
the channel walls. The method is developed for incompress- 
ible and linearized (y— — 1.0) compressible flow. The veloc- 
ity distribution along the channel walls is not arbitrary and 
the design method applies to elbows only. In reference 2 
a design method is developed for straight, symmetrical chan- 
nels with contracting or expanding walls. The method is 
developed for incompressible flow and the velocities are 
prescribed not as a function of arc length along the channel 
walls but as a function of circle angle in the transformed 
circle plane. A more general design is suggested in reference 3, 
but no attempt is made to develop and apply the method. 

In the present report a general method of design is devel- 
oped for two-dimensional, unbranched channels' with pre- 
scribed velocities as a function of arc length along the channel 
walls. The method is developed for both compressible and 
incoiiipressible, irrotational, honviscous flow and applies to 
the design of elbows, diffusers, nozzles, and so forth. Two 
types of compressible flow are considered: the general type 
with arbitrary value of y (1.4, for example) and the linearized 
type with y equal to —1.0. In general, if the prescribed 
velocity along one channel-wall differs from that along the 
other, the channel turns so that the downstream flow direc- 
tion is different from the upstream direction. This change 
in flow direction cannot be arbitrarily chosen but depends on 
the prescribed velocity distribution along the walls. Equa- / 
tions are developed for computing this change in flow direc- 
tion for an arbitrary prescribed velocity distribution with 
incompressible or linearized compressible flow. Two meth- 
ods of solution have been developed for the design method 
and are presented in separate parts of this report. In 
parti solutions are obtained by relaxation methods (ref. 4). 
This method of solution results in complete information 
concerning the distribution of flow conditions throughout the 

Velocity Distributions Along the Channel Walls. I — Relaxation Solutions” by John 
Velocity Distributions Along the Channel Walls. H— Solution by Green’s Function” by 


1 Supersedes NAOA TN 2593, "Design of Two-Dimensional Channels with Prescribed 
D. Stan Its, 1952, and NAOA TN 2595, "Design of Two-Dimensional Channels with Prescribed 
John D, Btanitx, 1052. 
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channel and, in addition, can be used to obtain nonlinear 
solutions for compressible flow with arbitrary values of y. 
In part II solutions are obtained by means of a Green’s 
function. This method of solution is limited to incom- 
pressible and linearized (7= — 1.0) compressible flow, but 
the method is more rapid than relaxation methods, provided 
information within the channel is not required. 

The design method reported herein was developed at the 
NACA Lewis laboratory during 1950 and is part of a doctoral 
thesis conducted with the advice of Professor Ascher H. 
Shapiro of the Massachusetts Institute of Technology. 

PART I 

GENERAL THEORY AND SOLUTION BY 
RELAXATION METHODS 

A general method of design is developed for two- 
dimensional, unbranched channels with prescribed velocities as 
functions of arc length along the channel walls. The method 
is developed for both incompressible and compressible, irro- 
tational, nonviscous flow. Two types of compressible flow 
are considered : the general type with arbitrary value for the 
ratio of specific heats y (1.4, for example), and the linearized 
type in which y is equal to —1.0. The solutions in part I 
of this report are obtained by relaxation methods and give 
complete information concerning the flow throughout the 
channel. Five numerical examples are given, including three 
elbow designs with the same prescribed velocity as a function 
of arc length along the channel walls hut with incompressible, 
linearized compressible, and compressible flow. 

THEORY OF DESIGN METHOD 

The design method is developed for two-dimensional chan- 
nels with prescribed velocities along the channel walls. The 
prescribed velocity is arbitrary except that stagnation points 
cannot be prescribed. This exception limits the design 
method to unbranched channels. 

PRELIMINARY CONSIDERATIONS 

0 

Assumptions. — The fluid is assumed to he nonviscous and 
either compressible or incompressible. The flow is assumed 
to be two dimensional and irrotational. 

The assumption of two-dimensional, nonviscous, irrota- 
tional motion limits the design method in practice to channels 
■with thin (negligible) boundary layers, such as exist near the 
entrance to the chhnnel or after a rapid acceleration of the 
flow through a contraction in the channel. Even if the 
^ boundary layer is thin, the design method is limited to (and 
finds its most useful application for) prescribed velocity 
distributions that, from boundary-layer theory, do not decel- 
erate fast enough to result in separation of the boundary 
layer, which separation alters the “effective” shape of the 
channel and completely changes the character of the flow. 

In some channels with fully developed turbulent boundary 
layers, the design method might be expected to yield results 
that are satisfactory, although approximate, because for this 


type of flow the rotational motion occurs primarily in regions 
dose to the channel walls. In ch arm el walls with thick or 
fully developed laminar boundary layers the design mothod 
cannot be used, because not only is the rotation of the flow 
important in most of the channel but, if the channel bends, 
important secondary flows develop that are not considered 
by the two-dimensional design method. 

Flow field. — The flow field of the two-dimensional chan- 
nd is considered to lie in the physical zy-plano where x and 
y are Cartesian coordinates expressed as ratios of a charac- 
teristic length equal to the constant channel width down- 
stream at infinity. (All symbols are defined in appendix A.) 

At each point in the channd (fig. 1) the velocity vector has 
a magnitude Q and a direction d where Q is the fluid velocity 
expressed as the ratio of a characteristic vdocity equal to the 
constant channd vdocity downstream at infinity. For con- 
venience, the vdocity Q is related to a vdocity g by 

i=Qs<t ( 1 ) 

where g is the velocity expressed as a ratio of the stagnation 
speed of sound and the subscript d refers to conditions down- 
stream at infinity. 

The flow direction 6 at each point in the channel is measured 
counterdockwise from the positive z-axis. From figure 1 

dx—ds cos 6 (2a) 

dy=ds sin 0 (2b) 

where ds is a differential distance in the direction of Q, that 
is, along a streamline. 



X 

Fiouei l— M agnitude and direction of valodty at point lnijr-plano. 
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Stream function and velocity potential. — If the condition 
of continuity is satisfied, a stream function p can be defined 
such that 


dp = pQ in 


(3) 


continuity and irrotational fluid motion expressed in terms of 
the transformed coordinates <p and p. 

Continuity. — The continuity equation expressed in terms 
of <p and p becomes (appendix B) : 


where p is the fluid density expressed as the ratio of a charac- 
teristic density equal to the stagnation density and where in 
is a differential distance measured normal to the direction 
of Q, that is, normal to a streamline. Along a streamline, 
in is zero so that from equation (3) the stream function p 
is constant. 

( 

If the condition of irrotational fluid motion is satisfied, a 
velocity potential <p can be defined such that 


1 /d log, p . d log, Q\ . b6 
p \ dip dp J'^P 


(5) 


Irrotational fluid motion. — The equation for irrotational 
fluid motion, expressed in terms of <p and p, becomes (ap- 
pendix B) : 

d log* Q b8 n /z>\ 

p — p d^~° (6) 


dcp—Q ds (4) 

Normal to a streamline, is is zero so that from equation (4) 
the velocity potential ip is constant. Thus lines of constant 
<p and p are orthogonal in the physical xy-plane. 

Outline of method. — Solutions of two-dimensional flow 
depend on known conditions imposed along the boundaries 
of the problem. In the inverse problem of channel design, 
the geometry of the channel walls in the physical zy-plane is 
unknown. This unknown geometry apparently precludes 
the possibility of solving the problem in the physical plane 
and necessitates the use of some new set of coordinates, that 
is, a transformed plane, in which to solve the problem. 
These new coordinates must be such that the geometric 
boundaries along which the velocities are prescribed are 
known in the transformed plane. It is also desirable, for 
mathematical simplicity, that the coordinate system in the 
transformed plane be orthogonal in the physical plane. A 
set of coordinates that satisfies these requirements is provided 
by <p and p, which are orthogonal in the physical zy-plane and 
for which the geometric boundaries are known constant 
values of p in the transformed y '/'-plane. . The distribution of 
velocity as a function of <p along these boundaries of constant 
p is known because, if 

Q=Q(*) 

is prescribed, equation (4) integrates to give 


Differential equation for distribution of velocity. — The 
second-order partial differential equation for the distribution 
of log, Q in the transformed y^-plane is obtained by differenti- 
ating equations (5) and (6) with respect to <p and p, respect- 


ively, and combining to eliminate 


b*9 

bipbp 


Thus, 


a 3 log, P a 3 log, q a log, p / a log, p a log, q \ 
a <(? a<p s h<p \ a^ a^> ) 

2 a log, q a log, P , 2 a 1 log, q n 
p ~a* a$ l_p W (7) 


Equation (7), together with a relation between p, Q, and g* 
determines the distribution of log, Q in the W'-plane for com- 
pressible flow with a given value of g d and for arbitrarily pre- 
scribed variations in log, Q along the boundaries of constant p. 

Density. — The density p is related to the velocity g by 
(ref. 5, p. 26, for example) 

<8a) 

which, from equation (1), becomes 

p =(i— (8b) 


<p=<p(s ) 

from which equations, 

Q=Q(v) 

The technique of the proposed method of channel design is 
therefore to obtain a differential equation for the distribution 
of velocity in the ^-plane. The velocity distribution ob- 
tained from the solution of this equation is then used to obtain 
the distribution of flow direction, from which distribution the 
channel walls in the physical ry-plane are obtained directly. 
The differential equation for the distribution of velocity in 
the v^-plano is nonlinear (for compressible flow with y other 
than —1.0) and is solved by numerical methods (relaxation 
methods). 

DIFFEItENTIAL EQUATION FOR DISTRIBUTION OF VELOCITY IN 
TRANSFORMED s^-PLANE 

The differential equation for the distribution of velocity in 
the transformed Wi-plane is obtained from the equations for 


Equation (8b) relates the density p to the velocity Q for a 
given value of g d . 

Incompressible flow. — For incompressible flow p is con- 
stant and equal to 1.0 so that equation (7) becomes 

b 3 log, Q b 2 log, Q /qn 

a.? 3 ■*" bp? u w 

Equation (9) determines the distribution of log, Q in the 
ipi/'-plane for incompressible flow. 

CHANNEL- WALL GEOMETRY 

After equation (7) or (9) has been solved to obtain the 
distribution of log, Q in the transformed ^-plane (for the 
arbitrary specified variations in log, Q with <p along the 
boundaries of constant p), the geometry of the channel 
walls in the physical zy-plane can be determined from the 
resulting distribution of flow direction 9. 
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Flow direction 0. — The distribution of flow direction 9 
along a stre amlin e (constant ip) is obtained from equa- 
tion (6), which integrates to give 

(10a) 

where the subscript ip indicates that the integration is taken 
along a line of- constant ip And where the constant of integra- 
tion is selected, to give a known value of 9 at one value of <p 
along each streamline. The integrand in equation (10a) is 
obtained from the distribution of- log, Q, which is .known 
from the solution of equation (7) or (9). 

The distribution of flow direction 9 along a velocity- 
potential line (constant <p) is obtained from equation '(5), 
which integrates to give 

(10b) 

where the subscript <p indicates that the integration is taken 
along a line of constant <p and where the constant of integra- 
tion is selected to give a known value of 6 at one value of 
\p along each velocity-potential line. As for equation (10a), 
the integrand in equation (10b) is known from the distribu- 
tion of log, Q obtained from the solution of equation (7) 
or (9). 

Channel-wall coordinates. — The variation in x along a 
line of constant ip in the p^-plane is given by 

bx_ /dx ds\ 
dp \ds dip/t 

which, combined with equations (2a) and (4), integrates to 
give 

x ~l^ dv (lla) 

Likewise, 

— JiW* (llb) 

(11c) 

. ' (lld) 

where the constants of integration are selected to give known 
values of a; or y at one value of p along each streamline or at 
one value of ip along each velocity-potential line. JEquations 
(11a) to (lid) determine the distribution of x and y in the 
transformed p '/'-plane or, which is the same thing, the shape 
of the strea mlin es and velocity-potential lines in the physical 
an/-plane. In particular, equations (11a) and (11c) when in- 
tegrated along the boundaries of constant ip in the pt/'-plane 
determine the shape of the channel walls. 

Turning angle. — In general, if the prescribed velocity dis- 
tribution along one channel wall differs from the distribution 
along the other wall, the channel deflects an amount A 9, 
which is the difference in flow direction far downstream and 
far upstream of the region in which the prescribed velocity 


distribution varies. In part II it is shown that for incom 
pressible flow the turning angle A 9 is given by 

A 9=9 d —9 u 



where the subscript u refers to conditions upstream at infinity 
-and where the subscripts 0 and 1.0 refer to the channe 
boundaries along which ip equals 0 and 1.0, respectively. A 
similar equation will be given later for the case of linearizec 
compressible flow. 

LINEARIZED COMPRESSIBLE FLOW 

The nonlinear differential equation (7) for the distributior 
of velocity in the p^-plane with compressible flow becomes 
linear and is considerably simplified if a linear variation ir 
pressure with specific volume (1/p) is assumed. This lineal 
relation between pressure and specific volume was first sug- 
gested by Chaplygin (ref. 6) in order to linearize the differ- 
ential equations for two-dimensional compressible flow ii 
the hodograph plane. 

Density. — If a linear variation in pressure with specific 
volume is assumed, the density p* is related to the velocity 
q* by (appendix C) 

P*= (l+2**)- 1/2 (13) 

where 

P*=£ lP (13a) 

and 

q*=faq (13b) 

where the constants ki and Jc 2 have been determined so tlmi 
values of p given by equation (13) equal the values of ^ 
given by equation (8a) for any two selected values of q (des- 
ignated by q a and q b ) . Thus, 



where p, and p» are determined by equation (8a) for the se- 
lected values of q a and g>, respectively. A discussion of the 
selection of q a and g> is given in appendix C. It will be 
noted that, if y is equal to —1.0, equation (8a) has the same 
form as equation (13). 

Stream function and velocity potential. — For the case o: 
linearized compressible flow it is convenient to define the 
stream function ip* and the velocity potential tp* by 


and 


dp*=p*q* dn 


(is; 

( 16 ) 


d<p*=q* ds 
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Continuity. — The continuity equation expressed in terms 
of <p* and yp* becomes (appendix D) 


Channel-wall coordinates. — The variation in x along a 
line of constant yp* in the £>*^*-plane is given by 


where 


or, conversely, 


b log« u | dd 


l+Vi+i* 5 



(17) 

(18) 
> (19) 


Irrotational fluid motion. — The equation for irrotational 
fluid motion, expressed in terms of <p* and yp*, becomes 
(appendix D) 


£>log„ n dd 

hyp* W 


( 20 ) 



which combined with equations (2a) and (16) integrates 
to give 


Likewise, 


r cos 0 , * 

x= — x- d<p 

Jr q 

(23a) 

C sin0 , ,* 
Jr p q 

(23b) 

C sin 0 , * 
y = —tt d<p 

Jr q 

(23c) 

f COS 0j,* 

v= „*„* 

Jr p q 

(23d) 


Differential equation for distribution of log,«. — The partial 
differential equation for the distribution of log, u in the 
p*^*_plane is obtained by differentiating equations (17) 
and (20) with respect to <p * and yp*, respectively, and com- 


bining to eliminate 


d 2 0 


Thus 


d 3 log, u , b 2 logs'U_ n 
+ 5^" =0 


( 21 ) 


Equation (21) determines the distribution of log, u in the 
y>*i/'*-plane for linearized compressible flow with a given 
value of gd and for arbitrarily prescribed variations in 
log, Q, related to log, u by equations (1), (13b), and (18), 
along the boundaries of constant yp*. Equation (21) is 
linear and is, like equation (9) for the case of incompressible 
flow, the equation of Laplace. Thus an incompressible 
flow- "solution for the distribution of log, Q in the y>^-plane is 
also a linearized compressible flow solution for the distribution 
of log, u in the p’^-plane. The transformation from the 
p^-plane is different, however, from the transformation from 
the <p‘V*-plane so that different channel shapes result in the 
sey-planc. r 

Flow direction 0 . — The distribution of flow direction 0 
along a streamline (constant yp*) is obtained from equa- 
tion (20), which integrates to give 

-S,*#** (22a > 


Likewise, the distribution of flow direction 0 along a velocity- 
potential line (constant <p*) is obtained from equa- 
tion (17), which integrates to give 

»— *■» 

Equations (22a) and (22b) for linearized compressible flow 
correspond to, and are used in the same manner as, equa- 
tions (10a) and (10b) for the usual type of compressible or 
incompressible flow. 


Equations (23a) to (23d) determine the distribution of x and 
y in the transformed p*^*-plane or, which is the same thing, 
the shape of the streamline and velocity-potential lines in the 
physical ay-plane. In particular, equations (23a) and (23c), 
• when integrated along the boundaries of constant yp* in the 
(pV*-plane, determine the shape of the channel walls. 
Equations (23a) to (23d) for linearized compressible flow 
correspond to, and are used in the same manner as, equa- 
tions (11a) to (lid) for the usual type of coihpressible or 
incompressible flow. 

Turning angle. — In part II it is shown that for linearized 
compressible flow the turning angle, or difference in flow 
direction far downstream and far upstream of the region in 
which the prescribed velocity distribution varies along the 
channel walls, is given by - 



(24) 


where A yp* is the value of yp* along the left boundary (channel 
wall) when faced in the direction of flow if the value of yp* 
along the right boundary is zero, and where the subscript 
A yp* refers to the boundary along which yp* is equal to A yp*. 


NUMERICAL PROCEDURE 

The channel design method in part I of this report was 
developed for three types of fluid flow: (1) compressible, 
(2) incompressible, and (3) linearized compressible. Al- 
though the numerical procedures of the design method are 
similar for each type of fluid, the procedures differ in detail 
and are therefore considered separately in this section. 

COMPRESSIBLE FLOW 

The numerical procedure for channel design with com- 
pressible flow (y=1.4, for example) is as follows: 

(1) The velocity is specified as a function of arc length 
along that portion of the channel walls over which the ve- 
locity varies 

’ 2 = 2 («) 
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< 2 = <300 (25) 

where s is arbitrarily equal to zero at that point along one 
channel wall where the velocity first begins to vary. 

(2) The channel-wall boundaries of the flow field in the 
transformed ^-plane are straight, parallel lines of constant 
4/ extending indefinitely far upstream and downstream 
between <p equals ± <= } where <p is arbitrarily equal to zero 
at that point on the channel wall at which s is equal to zero. 
The value of <p along the right channel wall when faced in the 
direction of flow (direction of positive <p ) is arbitrarily set 
equal to zero in which case the value of \p along the left 
channel wall (A\t) is obtained by integrating equation (3) 
across the channel at a position far downstream where flow 
conditions are uniform 

Aip—pd (26) 


(3) The distribution of log, Q as a function of <p along the 
boundaries in the p^-plane is obtained by integrating equa- 
tion (4) between limits so that 


?>= J** Qds =¥>(«) 


(27) 


which together with equation (25) gives the distribution of 
log, Q along the boundaries in the p^-plane 


log, Q=f (<p) (28) 

The integration indicated by equation (27) is carried out 
numerically for arbitrary distributions of Q as a function of s. 

(4) If the velocities prescribed along one channel wall 
differ from those along the other wall, the channel will, in 
general, turn the flow. This turning angle cannot be de- 
termined exactly for compressible flow until the channel 
design is completed. However, it will be shown that this 
turning angle is only slightly greater than that resulting for 
linearized compressible flow with the same prescribed 
velocity and with a suitable selection for q a and & in equa- 
tions (14a) and (14b). This latter turning angle for linearized 
compressible flow us given by equation (24), which can be 
integrated numerically for the arbitrary distribution of 
log, u=f (<p) corresponding to equation (28). 

(5) In order to solve equation (7) for the distribution of 
log, Q in the p^-plane, it is convenient to eliminate the 
density terms from equation (7) by means of equation (8b) . 
Thus, equation (7) becomes 


4$ imi +b51 ^ 1+40 (^j + 4D (*i2fcjy_ 0 

(29) 


where 


A=- 






'+1 

•-1 



Equation (29) is nonlinear, and it can be solved by relaxation 
methods (refs. 4 and 7, for example). A grid of equally 
spaced points, at each of which the value of log, Q is to be 
determined, is placed in the flow field between the channel- 
wall boundaries. The grid is extended upstream and down- 
stream sufficiently far- so that constant values of log, Q, are 
obtained across the channel by the relaxation methods. In 
the numerical examples to be presented six or eight grid 
spaces were used across the channel. In example HI the 
number of grid spaces was reduced from eight to four with 
negligible effect on the resulting channel design. The values 
of log, <2 at each grid point were relaxed to five significant 
figures. If the same velocity distribution is prescribed along 
both walls, the channel is symmetrical so that the velocity 
distribution in only one half of the channel need be deter- 
mined by relaxation methods. 

(6) After log, Q has been determined at each grid point in 
the ■pi/'-plano, the distribution of 0 is determined by equa- 
tions (10a) and (10b), which are integrated numerically. The 
constants of integration in equations (10a) and (10b) are 
determined to give a specified value of 0 at one point in the 
channel (far upstream, for example). The integrands in 
equations (10a) and (10b) are determined by numerical 
methods (tables I to VIE, ref. 4, for example) from the known 
values of p and log, Q at each of the grid points. If it is 
desired to know the flow direction along the channel-walls 
only, equation (10a) can be solved along the channel-wall 
boundaries ^=0 and \p=A\p only. If it is desired to know 
0 everywhere in the channel, the recommended procedure is to 
determine the variation in 0 along the mean streamline 
(4'= (Ayp)/2) by equation (10a) and to determine the 
variation in 0 along each velocity-potential line from the 
previously determined values on the mean streamline by 
equation (10b). 

(7) After the distributions of log, Q and 0 are known in the 
^-plane, the shapes of the streamlines and the velocity- 
potential lines in the physical zy-plano or, which is the same 
thing, the distributions of x and y in the transformed ip\p- 
plane are determined by the numerical integration of equa- 
tions (11a) to (lid). The constants of integration in these 
equations are determined so that specified values of x and y 
occur at one point in the flow field. The recommended 
procedure is to determine the variation in x and y along the 
mean streamline by equations (11a) and (11c) and to deter- 
mine the variation in r and y along each velocity-potential 
line for the previously determined values on the mean 
streamline by equations (lib) and (lid). If it is desired to 
know the x and y coordinates for the channel walls only, 
equations (11a) and (11c) can be solved along the channel- 
wall boundaries ^=0 and \p=A\p only. 


B= 1.0 
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INCOMPRESSIBLE FLOW 

The numerical procedure for channel design with incom- 
pressible flow (p=l) is s imil ar to that just outlined for 
compressible flow, but with the following differences: 

(1) The velocity is specified as a function of arc length by 
equation (25) alone. 

(2) The value of along the left channel wall (A \f/) is equal 
to 1.0 instead of the value given by equation (26). 

(3) The distribution of log, Q as a function of p along the 
channel-wall boundaries in the p^-plane is the same as that 
obtained from equations (25) and (27) and given by equa- 
tion (28) . 

(4) The turning angle AO of the channel is given by 
equation (12). 

(5) The distribution of log, Q in the p^-plane is obtained 
from the solution of equation (9) by relaxation methods. 

(6) After log, Q has been determined at each grid point 
between the channel-wall boundaries in the pi/'-plane, the 
distribution of 6 is determined by equations (10a) and (10b) 
as indicated previously for compressible flow, but with p 
equal to unity. 

(7) After the distributions of log, Q and 6 are known in the 
pi/'-plane, the shapes of the streamlines and velocity-potential 
lines in the physical zy-plane are determined by equa- 
tions (11a) to (lid) as indicated previously for compressible 
flow, but with p equal to unity. 

LINEARIZED COMPRESSIBLE FLOW 

The numerical procedure for channel design with line- 
arized compressible flow (7= — 1.0) is similar to that pre- 
viously outlined for compressible flow, but with the following 
differences: 

(1) The velocity g is specified as a function of arc length 
along the channel walls by g(s) or by g d and equation (25). 
For each prescribed velocity, there are an infinite number of 
linearized compressible flow solutions depending on the 
selected values of q a and g b in equations (14a) and (14b). 
However, for values of g a and g b within the range of g prescribed 
along the channel walls (and therefore everywhere in the 
channel), the solutions, that is, channel shapes, probably 
differ only in small detail. The best solution is that most 
nearly like the nonlinear compressible solution with arbitrary 
value of 7 (1.4, for example). In the numerical examples of 
this report it is shown that, if g a and g b are equal to the maxi- 
mum and minimum values of g, a good solution results, at 
least if the ratio of these prescribed velocities is not too large 
(2:1 in the numerical examples). On the other hand, if con- 
tinuity is to be satisfied for a gas with the correct value of 7 
(1.4, for example) upstream and downstream of the region 
of the channel in which the prescribed velocities vary, then 
g a and g b must equal g„ and g*. 

After ga and g b have been selected, the velocity distribution 
g(s) is expressed as g*(s) by equation (13b) where fa is given 
by equation (14b) so that 

g*=2 *(«) (30) 

The velocity g* is then expressed as u by equation (18) so 
that 

u=u{s) (31) 


In the particular case where the selected value of g n is equal 
to gi,, the value of fa is given by equation (C4b) in appendix O, 
where the significance of this particular case is also discussed. 

(2) The solution is obtained in the transformed p*^*- 
plane where p* and are defined by equations (16) and (15), 
respectively. If the value of </-* along the right channel wall 
when faced in the direction of g* is zero, the value of along 
the left wall (A^*) is obtained by integrating equation (15) 
across the channel at a position far downstream where flow 
conditions are uniform 

Af*=p d *q d * (32) 

(3) The distribution of log, u as a function of p* along the 
channel-wall boundaries in the p*^*-plane is obtained by 
integrating equation (16) between limits similar to those 
discussed previously for compressible flow so that 

P*= fV ds=p*(«) (33) 

which together with equation (31) determines the distribu- 
tion of log, u along the channel-wall boundaries in the p*^*- 
plane 

log, «=/(p*) (34) 

(4) The turning angle AO of the channel is given by equa- 
tion (24) . 

(5) The distribution of log, u in the p*^*-plane is obtained 
from the solution of equation (21) by relaxation methods. 

(6) After log, u has been determined at each grid point 
between the channel-wall boundaries m the p*</'*-plane, the 
distribution of 0 is determined by equations (22a) and (22b) 
in a manner similar to that outlined previously for compress- 
ible flow. 

(7) After the distributions of log, u and 0 are known in the 

p*i/'*-plane, the shapes of the streamlines and the velocity- 
potential lines in the physical ry-plane are determined by 
equations (23 a) to (23d) in a manner similar to that outlined 
previously for compressible flow. The velocities g* in equa- 
tions (23) are obtained from the known values of u, and the 
densities p* are given by equation (13) . , 

NUMERICAL EXAMPLES 

The channel design method has been applied in part I to 
the five examples listed below : 


Examples 

Type of channel 

Type of flow 

I 

Reducing section 

Incompressible 

n 

Converging section 

Incompressible 

m 

Elbow 

Incompressible 

, iv 

Elbow 

Linearized compressible 

V 

Elbow 

Compressible ( 7 — 1 .4) 


EXAMPLE I 

The first numerical example is the design of a reducing 
section' in a straight channel such that the upstream velocity 
is half the downstream velocity. The solution is for incom- 
pressible flow. 
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Fiqure 2. — Prescribed velocity distribution as function of arc length along channel wall for 
examples I, III, IV, and V. Equation (35). 

Prescribed velocity distribution. — The prescribed velocity 
as a function of arc length s along both ch ann el walls is 
given by 

Q= 0.5 (*<0) \ 

<3=|+~|^ (0<«<3.0)l (35) 

<3=1-0 (s<3.0) ) 

The prescribed velocity given by equation (35) is plotted in 
figure 2. 

Equation (35) together with equation (27) results in 
<p—0.5s . (*<0) \ 

0-5+R-R5 CO<*<3.oA (36) 

p= — 0.75+8 (s < 3.0) ) 


From equations (35) and (36), log, Q is a known function of «/>, 
which function is plotted in figure 3. 

Ee suits. — The results of example I are presented in 
figures 4 to 7. 

In figure 4, lines of constant velocity Q and flow direction 
6 are plotted in the transformed p^-plane. The flow direc- 
tion 6 is constant and equal to zero along the mean streamline 
(^=0.5), indicating that the center line of the channel is 
straight. The maximum absolute values of 9 occur along the 
channel walls. The solution is svmmetrical about the mean, 
streamline. The lines of constant Q and 8 are orthogonal. 


9 



Figure 3. — Prescrlbed^distrlbutionof log* Q as'fanction of *>[along channel walls for ox am pi o I, 
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♦ 



FiatntB Lines of constant x and t coordinates la transformed ^d-plane for example L Incompressible flow; prescribed velocity given In figure 2. 


X 



FiotJBE 6.— Streamlines and veloolty-potential lines on physical xp-plane tor example I. Incompressible flow; presodbed velocity given In figure 2. 


In figure 5, lines of constant a; and y are plotted on the 
transformed v^-plane. Along the mean streamline (^=0.5) 
the value of y is constant and equal to zero indicating, as 
before, that the center line of the channel is straight. The 
lines of constant x and y are orthogonal, and the system of 
curves forms a square network. The solution is symmet- 
rical. 


In figure 6, lines of constant <p and (velocity potential 
and streamlines, respectively) are plotted in the physical xy- 
plane. The shape of the channel walls is that required to 
result in the prescribed velocity distribution given by equa- 
tion (35) and plotted in figure 2. The downstream channel 
width is 1.0 by definition. The upstream channel width is 
2.0 in order that the upstream velocity be half the down- 


32100G — CD 12 
















162 


REPORT 1116 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

x 


4A rLE -s 1 ■ 'f i 'f . 2 1° ' 2 i 4 ' 2 ' S i - ^ -r 


\.z 


.8 


•4 


-.4 


-.8 


- 1.2 


n i i i r 


r i i i i i_ 



1.2 


.8 


0 y 


■A 


- 1.2 


i i l III I * I I I I L, J I 1 1 1 1 L 

1 1 ' "4 D L2 L6 2D 2^4 2D 3D 

x 


-L6 -1.2 -.8 -.4 0 

Fiodbe 7.— Lines of constant velocity Q and flow direction 0 in physical ay-plane for example L Incompressihle flow; prescribed velocity given in figure 2. 



Fioum; g, — Prescribed volocity distribution as function of arc kngth along channel wall for 
example IL Equation (38). 


stream velocity. As usual, the streamlines and velocity po- 
tential lines are orthogonal and, with equal increments of 
<p and \p, form a square network for incompressible flow. 

In figure 7, lines of constant Q and 8 are plotted in the 
physical ry-plane. The lines of constant Q and. 8 are orthog- 
onal. 

example n 

The second numerical example is the design of a converging 
section that f unn els the fluid from an infinite expanse into a 


straight channel of unit width. Far upstream the channel 
walls are straight and converge at a 90° angle. The solu- 
tion is for incompressible flow. 

Prescribed velocity distribution. — The prescribed velocity 
as a function of arc length s along both channel walls is given 

by 


«-^5> ( ^ 0) . 

r&- O' (0 ^ ‘ 4> P> 

Q= 1.0 (s>4) ) 


The prescribed velocity given by equation (37) is plotted in 
figure 8. 

Equation (37) together with equation (27) results in 


r k *-( 1 -f) (>so > 

- 1 1 £! 1 (1 3\ 

^ v S ' 2v 2 8 \2 tt 2/ 3 

n(H? <0Ss£4) 

¥>=^—2+8 (s>4) 


(38) 


From equations (37) and (38), log, Q is a known function of 
•<p, which function is plotted in figure 9. 
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Figure 9.— Prescribed distribution of log, Q as function of along channel walls for example n. 

I 

Results. — The results of example II are presented in 
figures 10 to 12. 

In figure 10, lines of constant velocity Q and flow direction 
0 are plotted in the transformed ■pi/'-plane. The flow direc- 
tion 0 is constant and equal to zero along the mean stream- 
line (i/'=0.5), indicating that the center line of the channel 
is straight. The solution is symmetrical about the mean 
streamline. As for example I, the lines of constant Q and 0 
are orthogonal. 

In figure 11, lines of constant <p and \f/ are plotted in the 
physical zy-plane. The shape of the channel walls is that 


required to result in the prescribed velocity distribution 
given by equation (37) and plotted in figure 8. As usual, the 
streamlines and velocity-potential lines are orthogonal and, 
for incompressible flow with equal increments of <p and V', 
form a square network. 

In figure 12, lines of constant Q and 0 are plotted in the 
physical zy-plane. The lines of constant Q and 0 are orthog- 
onal. 

example m 

The third numerical example is the design of an elbow for 
which the upstream velocity is half the downstream velocity. 
The prescribed velocities are such that no deceleration occurs 
anywhere along the channel walls. The solution is for in- 
compressible flow. 

Prescribed velocity distribution. — Along both walls up- 
stream of the elbow the velocity Q is equal to 0.5, and along 
both walls downstream of the elbow Q is equal to 1.0. The 
transition from Q equals 0.5 to 1.0 along both walls of the 
elbow will be the prescribed velocity distribution as a func- 
tion of arc length given by equation (35) for example I and 
plotted in figure 2. In terms of log, Q as a function of <p, this 
prescribed velocity distribution is given by equation (36) 
and is plotted in figure 3. Although this velocity distribu- 
tion is the same for both walls, the distribution on the outer 
wall (wall with larger radii of curvature) is shifted in the 
positive cp direction an amount equal to 2.25 relative to the 
distribution on the inner wall. Thus, a velocity difference 
exists on the two walls at equal values of <p, as shown in 
figure 13. The greater this difference in velocity and the 
greater the range in <p over which velocity differences exist, 
the greater is the elbow turning angle. For the prescribed 
velocity distribution given in figure 13, the elbow turning 
angle given by equation (12) was 89.37° compared with a 
value of 89.36° obtained from the relaxation solution. 

Results. — The results of example HI are presented in 
figures 14 to 16 and in tables I and II. (The numerical 



1.0 

.8 

.6 

.4 

.2 

.0 


Figure 10.— Lines of constant velocity Q and flow direction 8 In transformed p^-plane for example IL Incompressible flow; prescribed velocity given In figure 8. 
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Fiqtoe 11. — Streamlines and velocity-potential lines in physical xy-plane for example IL Incompressible flow; prescribed velocity given in figure 8. 
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Figure 13.— Prescribed distribution ol log. Q as function of along channel walls for examples 

HI, IV, and V. 


results for examples HI, IV, and Y are tabulated in tables 
I to VI to enable a detailed comparison of the three elbow 
designs with the same prescribed velocity Q distribution as 
a function of arc length but with incompressible (example 
111), linearized compressible (example IV), and compressible 
(example V) flow.) 

_ In figure 14, lines of constant Q and 6 are plotted in the 
^i/'-plane. The flow direction 8 varies along the mean 
streamline (^=0.5), -indicating that the channel is curved. 
The solution is unsymmetrical. As for examples I and II, 
the lines of constant Q and 6 are orthogonal. 

In figure 15, lines of constant <p and \p are plotted in the 
physical xy-plane. The shape of the channel walls is that 
required to result in the prescribed velocity distribution 
given by equations (35) and (36) and plotted in figures 2 
and 13. The upstream channel width is twice the down- 
stream width in order that the upstream velocity be half 
the downstream velocity. It is interesting to note that, 
before curving in the direction of the elbow turning angle, 
the inner wall first curves in the opposite direction. This 
behavior of the inner-wall geometry is necessary in order to 


maintain the prescribed constant velocity along the outer 
wall where the velocity would otherwise decelerate because 
of the necessary curvature in the direction of elbow turn- 
ing. This feature of the elbow geometry will also bo noted in 
examples IV and V. As usual, the streamlines and velocity- 
potential lines are orthogonal and, for equal increments of 
<p and \p, form a square network. 

In figure 16, lines of constant Q and 8 are plotted in the 
physical sy-plane. The lines of constant Q and 0 are 
orthogonal. 

EXAMPLE IV 

The fourth numerical example is the design of an elbow 
with the same prescribed velocity Q, as a function of arc 
length, used in example HI but for linearized compressible 
flow ( 7 = — 1 . 0 ). 

Prescribed velocity distribution. — The prescribed velocity 
distribution Q is the same as that for example IH and with 
gi equal to 0.80176. The variation in Q with s along one 
channel wall is plotted in figure 2 . The values of g„ and g» 
in equations (14a) and (14b) are equal to g« and g d , or 
0.40088 and 0.80176, respectively. For these values of g„ 
and gi and for the prescribed velocity distribution with 
linearized compressible flow, the elbow turning angle given 
by equation (24) was 104.08° compared with a value of 
104.07° obtained from the relaxation solution and a value of 
89.36° obtained for incompressible flow (example HI). 

Besults. — The results of example IV are presented in 
figures 17 to 19 and in tables HI and IV. 

In figure 17, lines of constant g and 8 are plotted in the 
transformed ¥>*^*-plane. The solution is unsymmetrical 
and the lines of constant g and 8 are orthogonal. 

In figure 18, lines of constant <p*/A\p* and rp*/A\p* are 
plotted in the physical xy-plane (where the constant Arp* 

, is given by equation (32) and is equal to 0.73782 for q d 
equal to 0.80176). The shape of the channel walls is that 
required to result in the prescribed velocity distribution used 
in example HI but with linearized compressible How and 
for gi equal to 0.80176. From continuity considerations 
the upstream channel width is 1.5385 times the downstream 
width. As in example IH, the inner wall of the elbow first 



Figure 14. — Lines of constant velocity Q. and flow direction 6 In transformed p^plane for example HL Incompressible flow; prescribed velocity given in figures 2 and 13. 



DESIGN OP TWO-DIMENSIONAL CHANNELS WITH PRESCRIBED VELOCITY DISTRIBUTIONS ALONG CHANNEL WALLS 167 


1 . 0 - 
,5- 
0 - 
-.5- 
-1.0- 
y -1.5 - 
- 2 . 0 - 
-2.5- 
-3.0- 
-3.5- 
-4.0- 


X 

-1.0 -.5 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 . 

— j i 1 1 1 nr =1 1 1 1 1 


* 



l I ! I 1 1 1 I I I 1 

-1.0 -.5 0 .5 1.0 1.5 2.0 2.5 *3.0 3.5 4.0 

x 

Fiqtjee 15.— Streamlines and velocity-potential lines In physical ip-plane for example KL Incompressible flow; pre*cribed_veloolty given In figures 2 and 13. 
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Figure 16. — Lines of constant velocity Q and flow direction 6 in physical xy-plane for example HL Incompressible flow; prescribed velocity given In figures 2 and 13. 
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ip'/Aty* 

FiauBE 17.— -Lines of constant velocity q and flow direction Q in transformed V*-plane for example IV. Linearlxed compressible flow; prescribed velocity as function of are length along channel 

walls same as for example m (fig. 2) and with q* equal to 0.S0176. 


turns in the opposite direction to the elbow turning angle. 
As usual, the streamlines and velocity-potential lines are 
orthogonal. 

In figure 19, lines of constant g and 0 are plotted in the 
physical ay-plane. The lines of constant g and 0 are not, in 
general, orthogonal. 

example v 

The fifth numerical example is the design of an elbow with 
the same prescribed velocity Q, as a function of arc length, 
used in examples HI and IV but for compressible flow 
(7=1.4). 

Prescribed velocity distribution. — The prescribed velocity 
distribution <2 is the same as that for examples III and IV 
but with (ji equal to 0.79927. The variation in Q with s 
along one channel wall is plotted in figure 2 . 

Besults. — The results of example V are presented in 
figures 20 and 21 and in tables V and VI. 

In figure 20, lines of constant <p/Ai{/ and are plotted in 
the physical ay-plane (where the constant Af is given by 
equation (26) and is equal to 0.71054 for g 4 equal to 0.79927). 
The shape of the channel walls is that required to result in 
the prescribed velocity distribution used in examples m and 
IV but with compressible flow ( 7 = 1 . 4 ) and for g«j equal to 
0.79927. The upstream channel width is 1.5412 times the 
downstream width, and the turning angle is 105.31° com- 
pared with 104.07° for linearized compressible flow (example 
IV) and 89.36° for incompressible flow (example HI). The 
streamlines and velocity-potential lines are orthogonal. 

The shape of the elbow for compressible flow (example V, 
fig. 20 ) is nearly the same as the shape of the elbow for 
linearized compressible flow (example IV, fig. 18). There- 
fore, in figure 21 the contours of the walls for both examples 
are compared. The difference in contours is very small and 
it is concluded that, if a nonviscous gas with arbitrary 7 (1.4, 
for example) were to flow through a channel designed for 
linearized compressible flow ( 7 = — 1 . 0 ) , the resulting velocity 
distribution along the channel walls would be nearly the 
velocity distribution prescribed for the linearized compress- 
ible flow, at least if the linearized flow were selected (by 
the choice of g a and g & ) so that the densities were equal for 
both types of flow at the maximum and minimum velocities 


and if the ratio of these prescribed velocities is not too large 
(2 : 1 in the numerical examples) . This conclusion is impor- 
tant because the design method for linearized compressible 
flow is considerably faster than the design method for com- 
pressible flow with 7 other than — 1 . 0 . 

PART n 

SOLUTION BY GREEN’S FUNCTION 

In part II a method of solution for the design of two- 
dimensional channels with prescribed velocity distributions 
along the walls is developed by means of the appropriate 
Green’s function. The method applies to incompressible and 
linearized compressible, irrotational flow. One numerical 
example is presented for an accelerating elbow with linear- 
ized compressible flow and with the same prescribed condi- 
tions as example IV of part I. 

METHOD OF SOLUTION 

The method of solution by Green’s function is in conjunc- 
tion with a formula derived from Green’s theorem. 

PRELIMINARY CONSIDERATIONS 

Stream function \F. — In part II it is convenient to define 
the stream function by tF, where for incompressible flow 

df (39a) 


and for linearized compressible flow ( 7 =— 1 . 0 ) 


d^r— 



(39b) 


For both types of flow tF varies from zero along the right side 
of the channel, when faced in the -direction of flow, to ir /2 
along the left side. 

Velocity potential 4>, — In part II it is convenient to define 
the velocity potential by 4>, where for incompressible flow 


d^=~ dip 


(40a) 
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Figure 19. — Lines of constant velocity q and flow direction 6 In physical xy-plane for example IV. Linearised compressible flow; prescribed velocity as function of arc length along channel 

walls same as for example HE (fig. 2) and with q i equal to 0.80176. 
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Fwote 20.— Streamlines and velodty-potantlal Unas In physical ry- plane for eiampla V. Compressible flow (yl.4); prescribed velocity as function of aro length along ohannol walls same as 
a for eiamples HI and IV (fig. 2) but with gj equal to 0.79927. 
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Figure 21.— Comparison of channel-wall shapes for compressible flow (example V) with y equal to 1.4 and for linearised compressible flow (example TV) for same prescribed velocity as function 

of arc length along channel walla (fig. 2). 
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and for linearized compressible flow 

a 2A^ 


(40b) 


Channel-wall coordinates. — From part I the distribution 
of channel-wall coordinates x and y along the boundaries of 
constant 4 1 ' equal to 0 and ir/2 in the transformed i^-plane is 
given by 

2 ,* f cos ‘ 

x=—A\p* 

TT J •i 


and 


1* q 


d<f> 


(41a) 


2 ,* C sm 0 

y=— — ■*— d$ 

■x Jv q 


(41b) 


for linearized compressible flow, and for incompressible flow 
is given by 

cos 9 


and 


-2f 

irjn 


2 f si 

y=~ — 
vj-t ' 


<t Q 

sin 9 








(42a) 


(42b) 


where the constants of integration are selected to give known 
(specified) values of x or y at one value, of $ along each bound- 
ary. Because g* and Q are known functions of $ from the 
prescribed velocity as a function of arc length along the 
channel walls, the shape of the channel walls in the physical 
ay-plane is given by equation (41) or (42) if 9 is determined 
as a function of $ along the channel walls. In part II the 
solution for 9 as a function of $ along the channel walls in the 
$4 f -plane is obtained by Green’s function. 

SOLUTION BY GREEN'S FUNCTION 

Continuity. — From part I the continuity equation becomes 
in the transformed $fp-plane 


d log, V . d9 

where for incompressible flow 

V=Q 

and for linearized compressible flow 


=0 


1+Vl + y* 


(43a) 


(43b) • 


(43c) 


Irrotational motion. — From part I the equation for irro- 
tational motion becomes in the transformed ^-plane 


d log, V d9 






=0 


(44) 


Integral equation for 0($ O ,M' O ). — From equations (43a) and 
(44) 


V9 , d*e 


(45) 


so that from appendix E the value of 9 at a point ('f> 0) 4 r < ,) 


within, or on, the channel walls in the transformed 'I'l'-plano 
is given by the integral equation 

*<*-.**>-i 7 K e (46) 


where the subscripts 0 and — refer to the channel-wall bound- 

A 

aries along which is 0 and respectively, and 0 is the 
Green’s function of the second kind for the channel, which 
is an infinite strip of width ^ extending in the ^-direction to 
± 00 . 

Green’s function Q.—j The Green’s function of the second 
kind 6 for the infinite channel in the -plane is given along 

the channel-wall boundaries (4r equals 0 and ^ ) by (appen- 
dix F) 

GW 1 = —log, [cosh 2 ($—$„) — cos 2 (^ — 'I',)] (47) 

where ($, t?) is any point on the channel-wall boundary and 
($o,'3! r «) is the point in the channel or on the boundary Jat 
which 9 is to be determined. 

Numerical integration for 0($„,4 r „). — From equations (46) 
and (47) 

2v0($ o ,^„)= J ^|^|^log. [cosh 2 ($-$,)— 


sin 2 '®', 


,]} r d ($-$,)- 


/:{ 


log, [cosh 2 ($—$„) - 


cos 2 


j}od ($-$,) 


(48) 


in which the independent variable of integration has. been 
changed from to d($— $„) so that the origin, for purposes 
of integration, lies at 4> 0 rather than $=0. If for small 


changes in ($—$,), that is, for small A<J>, the term 


aiog.y 

d4> 


may be considered constant and equal to its average value 
over the interval A4>, then 

d log,y _ Alog,F 
5$ A4> 

and equation (48) becomes 

» (A lov VC (*-*.)+A* 

log, [cosh’ 


A$ 

sin 2 ’S' 0 ]'d(<l>— 4>, 


i- 


- (a \oz.V C (*-•*■ J+a* 

log ' [“**(*-*•)■ 


cos 2 Sk 


,] d($— 4> 0 )j° 


(49) 
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where the summation sign is understood to mean that the 
quantity within the braces is summed over the entire range of 
(<£—$<,) between ± <*> . 

Equation (49) determines 6 at any point in the flow field 
(channel). For a point ($„,'& 0 ) on the channel walls ik, is 
equal to 0 or v/2 and the integrands in equation (49) become 

2 log, cosh |(4>— 4>„)| 
or 

2 log, sinh |($— $,)| 


so that equation (49) becomes 
'A log,F 


J?~. K^t 1 "X-CtS 2 ")] (50a) 


where 


AJ= J(4_4,.)+A» — I&-Q). 


(50b) 


(2) Compute V as a function of s from equations (43b) and 
(51a) for incompressible flow or from equations (13b), (14b), 
(43c), and (51b) for linearized compressible flow 

V=V(s) ' (52) 

(3) Compute 4> as a function of s from equations (4) and 
(40a) for incompressible flow or from equations (16), (32), 
(40b), and (51b) for linearized compressible flow. In equa- 
tion (32) p d * is obtained from equations (8a), (13a), and (14a). 
For arbitrary distributions of Q or g equation (40a) or (40b) 
is integrated numerically by using, for example, Simpson’s 
one-third rule. Thus 

4>=4>(s) (53) 

(4) From equations (52) and (53) V and 4> are known 
functions of s so that 

V=F($) (54) 


Il=cc if 4V=0 1 
2 

* *‘=i 

► 

7„=aif ¥„=! 

!„=/ Sif\k o =0 . 

where 

a= ± log, cosh |($— $ 0 )| d |($— $„)| 

f l(+— M 

j9= ± J fl log, sinh |($— 4*„)| d |($— $„)| 


(50c) 


(50d) 

(50e) 


where the + signs apply for positive values of ($—$„) and 
the — signs apply for negative values of ($—$„). Methods 
of evaluating a and /3 are given in appendix G, and tabulated 
values are given for a wide range of |($— $„)| in table V1L. 
Equation (50a) determines 0($„4 r <1 ) at any point on the 
channel-wall boundaries. Thus from equations (41a) and 
(41b) or (42a) and (42b) the coordinates for the channel-wall 
shape in the physical zy-plane can be determined. 


NUMERICAL PROCEDURE 

The numerical procedure for the channel design solution 
by Green’s function is the same, except for minor details, 
for incompressible and linearized compressible flow. The 
stepwise procedure is outlined as follows: 

(1) For incompressible flow the velocity Q and for linear- 
ized compressible flow the velocity g, or which is the same 
thing the velocity Q and the constant downstream velocity 
g i, are specified as functions of arc length along the channel 
walls 

Q=Q(s) ■ (51a) 

or 

g=g(s) (51b) 

where s is arbitrarily equal to 0 at that point along one chan- 
nel wall where the velocity first begins to vary. 


Thus V is a known function of 4> along the channel-w'all 
boundaries in the transformed 4>4'-plane. 

(5) If the prescribed velocity distribution along one wall 
is different from that along the other, the channel will, in 
general, turn the flow. This turning angle AO is given by 
equation (H5) in appendix H. If the turning angle is unsat- 
isfactory, a new distribution of velocity as a function of s 
(eqs. (51a) and (51b)) is prescribed and steps (1) to (5) 
repeated until the desired value of A0 is obtained. Equation 
(H5) is integrated numerically by using Simpson’s one-third 
rule, for example, and equation (54) . 

(6) The channel-wall boundaries are straight parallel 
lines of constant 4' equal to 0 and ir/2, and extending to ± » 
in the ^-direction. Along these boundaries of constant 4', a 
series of equally spaced points are located at each of which 
the flow direction 0 and the r,y-coordinates of the channel 
walls will be determined by numerical integration. In 
order to use the tables of a and jS presented in this report, 
the point spacing A<3? must be an even multiple of tt/ 24. Thus 
the smallest point spacing -r/24 is equal to K» of the channel 
width ( 71 -/ 2 ). For a particular prescribed velocity distribu- 
tion along the channel walls the accuracy of the solution 
increases, and so does the amount of computing, as the 
point spacing is reduced. The error for a given point spacing 
depends on the prescribed velocity distribution, and its order 
of magnitude is given by the leading term of the error series 
of the formula used for numerical integration (table \'lll, 
ref. 4, for example). For the numerical example presented 
in part IE of this report the point spacing A$ was n-/12. 
From equation (54) 


A log, F~ _ (log, y)s+A$— (log, V)i 
■ A$ A$ 


(55) 


where the subscripts $ and 4>-bA4> refer to adjacent points 
along the channel boundaries. 

(7) The value of 6 at each point (^ 0 ,^o) on the channel- 
wall boundaries is obtained from equation (50a) in which 
(A log, F)/A$ is given by equation (55) and A I is given by 
equations (50b), (50c), and table V1L. Note that in equation 
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(50a) the origin has been moved to by changing from $ 
to ($— $ 0 ) . Thus the value of (A log, V) /A3? for a given value 
of ($—$„) varies with $>„. 

(8) The physical ^-coordinates at each point on the 
channel-wall boundaries are obtained by the numerical 
integration of equations (42a) and (42b) for incompressible 
flow, or equations (41a) and (41b) for linearized compressible 
flow where A\p* is given by equation (32). The constants of 
integration in equations (41) and (42) are selected to give 
known values of x and y at upstream or downstream positions 
where flow conditions can be considered uniform. 

NUMERICAL EXAMPLE 

The channel design method of part II has been applied to 
the design of an elbow for the same conditions as example 
IV of part I. The design is for an accelerating elbow with 
no local decelerations of the prescribed velocities along the 
channel walls and with linearized compressible flow. 

'Prescribed velocity distribution. — The prescribed velocity 
distribution along the channel walls is the same as that for 
example TV of part I. The prescribed velocity as a function 
of $ is plotted in figure 22. 

Results. — As indicated in table V1LL, the elbow design 
resulting from the prescribed velocities given in figure 22 
is the same as that obtained by relaxation methods (fig. 21) 
for the same prescribed conditions (example IV, part I). 

The solution obtained by Green’s function (part II) 
required one experienced computer 3 days, whereas the 
solution by relaxation methods (part I) required about 
10 days. The relaxation solutions provide additional 
information, such as the distribution of velocity across the 
channel; but for the most part this additional information 
is of secondary importance, and the design of channels by 
Green’s function is more rapid and therefore to be preferred 
over the design by relaxation methods. 


4 



F joube 22. — Variation In prescribed values of log. V with 4- along channel walls of numerical 

example in pert IL 


SUMMARY OF RESULTS AND CONCLUSIONS 

A general method of design is developed for two- 
dimensional unbranched channels with prescribed velocities 
as a function of arc length along the channel walls. The 
method is developed for both compressible and incompress- 
ible, irrotational, nonviscous flow and applies to the design 
of elbows, diffusers, nozzles, and so forth. Two types 
of compressible flow are considered : the general type with 
arbitrary value for the ratio of specific heats y (1.4, for 
example) and the linearized type in which y is equal to 
— 1.0. In part I solutions are obtained by' relaxation 
methods on a transformed plane the coordinates of which 
are the streamlines and velocity-potential lines in the physi- 
cal plane; in part II solutions are obtained by a Green’s 
function. The method of solution in part I gives complete 
information concerning the flow throughout the channel, 
whereas the method of solution in part II gives the channel- 
wall coordinates only. 

Five numerical examples are given in part I and the results 
are presented by (1) lines of constant velocity and flow 
direction or lines of constant physical coordinates in the 
transformed plane and (2) streamlines and velocity-potential 
lines or lines of constant velocity and flow direction in the 
physical plane. Among the fiye -examples are three elbow 
designs for the same prescribed velocity as a function of 
arc length along the channel walls but with incompressible, 
linearized compressible, and compressible flow. The numer- 
ical results of these three elbow designs are tabulated to 
enable a detailed comparison of the three designs. 

The shapes of the elbows for compressible flow and for 
linearized compressible flow are very nearly the same; and it 
is concluded that, if a- nonviscous gas with, arbitrary y 
(1.4, for example) were to flow through a channel designed 
for linearized compressible flow (y= — 1.0), the resulting 
velocity distribution along the channel walls would be 
nearly the. velocity distribution prescribed for the linearized 
compressible flow. This conclusion is important because 
the design method for linearized compressible flow is con- 
siderably faster than that for compressible flow. 

One numerical example is presented in part II for an 
accelerating elbow with linearized compressible flow. The 
elbow shape obtained from the solution by Green’s function 
in part II is the same as that obtained from a solution by 
relaxation methods in part I for the same prescribed con- 
ditions. The time, required for the calculations was con- 
siderably less for the solution by Green’s function. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 25, 1951 
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APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 


A, B, C, D 
A, B 

Bi, Bt, . . . 

c 

Q 

I 

h 

h 

l 

n 


V 

Q 


i 

q* 


r 


s 


u 

V 

W, Wi, Vh 


V 


z 

g 

a 

p 

7 

A 

S 


coefficients, equation (29) 
arbitrary constants, equation (Cl a) 
Bernoulli's numbers 
constant, equation (E3) 

Green’s function of the second kind, equa- 
tions (E2) and (47) 
integral (a or |S) 
coefficient, equation (14a) 
coefficient, equation (14b) 
length of closed boundary 
distance in zy-plane measured normal to 
direction of flow (expressed as ratio of 
characteristic length equal to channel 
•width downstream at infinity) 
static pressure (expressed as ratio of stag- 
nation density multiplied by stagnation 
speed of sound squared) 
velocity (expressed as ratio of characteristic 
velocity equal to constant channel veloc- 
ity downstream at infinity) 
velocity (expressed as ratio of stagnation 
speed of sound) 

velocity used in linearized compressible 
flow and related to q by equation (13b) 
distance from any point in ^ik-plane to 
point ($ 0 , T„) at which logarithmic singu- 
larity exists 

distance in xy-plane measured along direc- 
tion of flow (expressed as ratio of char- 
acteristic length equal to channel width 
downstream at infinity) 
velocity parameter related to q* by equa- 
tion (18) 

velocity parameter defined by equations 
(43b) and (43 c) for incompressible and 
linearized compressible flow, respectively 
complex functions defined by equations 
(F3), (Fla), and (F2a), respectively 
Cartesian coordinates in physical plane 
(expressed as ratios of characteristic 
length equal to channel width down- 
stream at infinity) 
complex coordinate, equation (Fib) 
conjugate of z 
integral, equation (50d) 
integral, equation (50e) 
ratio of specific heats 
finite increment 
increment of 


6 

Ad 

P 

P * 

$ 


<P and <p * 


¥ 


ip and ip* 
Ap* 


a 

Subscripts : 
a,b 


d 

o 

u 

Alp* 


($-$„) 


($-$ 0 )+A$ 


0 

1.0 

7 T 
2 


flow direction in physical cry-plane (meas- 
ured in counterclockwise direction from 
positive a- axis) 

channel turning angle, equation (12) 
density (expressed as ratio of stagnation 
density) 

density in linearized compressible flow and 
related to p by equation (13a) 
velocity potential used as Cartesian co- 
ordinate in transformed 'fik-plane and 
related to <p or <p* by equation (40a) or 
(40b), respectively 

velocity potential for incompressible and 
linearized compressible flow, respectively, 
equations (4) and (16) 
stream function used as Cartesian coordi- 
nate in transformed $\k-plane and related 
to \p or ip* by equation (39a) or (39b), 
respectively 

stream function for incompressible and 
linearized compressible flow, respectively, 
equations (3) and (15) 

boundary value of ip*, for linearized com- ' 
pressible flow, along left channel wall 
when faced in the direction of flow, 
equation (32) 

any harmonic function in ^T-plane 


quantities related to two velocities (q a and 
q t , respectively) for which density given 
by equation (8a) is equal to density p 
given by equations (13), (13a), and (13b) 
conditions downstream at infinity 
point in ^T-plane at which d is determined 
conditions upstream at infinity 
left channel wall, when faced in direction of 
flow, along which ip* is equal to Aip* 
point at ($—$„) on either channel-wall 
boundary 

point at [($—$„) +A3>] on either channel- 
wall boundary 

along lines of constant <p, ip, <p*, and ip*, re- 
spectively 

right channel wall, when faced in direction 
of flow, along which T, ip, or ip* is equal 
to 0 

left channel wall, when faced in direction 
of flow, along which \ p is equal to 1.0 
left channel wall, when faced in direction 

of flow, along which ik is equal to ~ 
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APPENDIX B 

EQUATIONS OF CONTINUITY AND IRROTATIONAL FLUID MOTION IN TERMS OF TRANSFORMED <p, ^-COORDINATES 


Consider the two-dimensional irrotational motion of a 
fluid particle in the physical ry-plane. The fluid particle is 
defined by adjacent streamlines (constant p) and velocity- 
potential lines (constant <p) spaced Sn and 5a apart as 
indicated in figure 23. The velocity Q is parallel to the 
streamlines and normal to the velocity-potential lines. 

Continuity. — From continuity considerations of the fluid 
particle in figure 23 


^(pQ5n)=0 


or 


a log, P , b log, Q 


a 8 


a 8 


__L d(to) _ Q 
Sn a a 


(Bl) 



Figure 23. — Fluid par tide bounded by streamline* and velocity-potential lines in physical 

rp-plane. 


But, from geometrical considerations (ref. 6, p. 167, for 
example) 


and 


Sn 

b(Sn) 

bs 

_be 

bn 

(B2a) 

JL 2>(5a)_ 
5a bn 

_be 

bs 

(B2b) 


so that equation (Bl) becomes 

a log, P , a log, Q 


be 


a s 


bs 


bn 


or 


a log, p dtp .b log, Q dtp . bO d^ 

b)tp 


=0 


ds ' btp ds ' bp dn 
which, combined with equations (3) and (4), becomes 


i / a log, p a log, q \ ae 
pV be "** bv J^bp 


(5) 


Equation (6) is the continuity equation expressed in terms 
of ^-coordinates. ' 

Irrotational fluid motion. — For irrotational motion of the 
fluid particle in- figure 23 • 


or 


£(«*>“ 0 

a log, Q 1 £>(5a) 
bn ' Ss bn 

But, from equations (B2b) and (B3) 

a log, Q be_ 


= 0 


(B3) 


an 


a« 


o 


or 


a log, Q dp_b6_ dtp 
bp dn btp ds 

which, combined with equations (3) and (4), becomes 

a log, Q be _ n 

p bp btp 0 

Equation (6) is the equation for irrotational fluid motion 
expressed in terms of the ^^-coordinates. 


( 6 ) 


APPENDIX C 

RELATION BETWEEN VELOCITY AND DENSITYIASSUMING LINEAR VARIATION IN PRESSURE WITH SPECIFIC VOLUME 


The approximate, linear relation between pressure p and 
specific volume 1/p first suggested by Chaplygin (ref. 6) is 
given by 

p=A-~ (Cla) 

p 

from which 

t=7' (01b) 


where A and B are arbitrary constants. 

If p denotes the static pressure expressed as a ratio of the 
stagnation density multiplied by the stagnation speed of 
sound squared, Bernoulli’s equation is 


— +qdq=0 
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which combined with equation (Clb) integrates to give the 
approximate relation between velocity and density 

B a 1 

2 ^ 5 — ■ 7 j-= constant (C 2 ) 

< 

For convenience equation (C 2 ) can he written as 



_L_0*s— 1 



P 


or 

P*=(l+2**) _l/J 

(13) 

where 

p*=kip 

(13a) 

and 

q*=kzq 

(13b) 


The constants ki and kz replace the two arbitrary constants 
in equation ( 02 ), and their values are determined so that 
for any two arbitrary values of q (designated by q a and q b ) 
the values of p given by equation (13) equal the values of p 
given by equation ( 8 a). Thus the values of p given by 
equation (13) for q equal to q a or q b are correct; for all other 
values of q the values of p are approximate. The constants 
ki and k 3 are determined from the conditions 


Pa* — kip a 

q a *=*kzq a 

Pt,*—kip b 

qt,*=kzq b 


(C3) 


and 



(14b) 


where p a and p b are determined by equation ( 8 a) for the 
selected values of q a and q b , respectively. 

The values of q a and q b might, for example, be selected to 
equal the maximum and minimum values of q (which values 
of q must occur on the channel walls and are therefore known) . 
Also, the values of q a and q b might be selected to equal the 
upstream and downstream velocities q u and q d . In this case 
the upstream and downstream channel widths would then 
satisfy continuity for a gas with the correct value of y (1-4, 
for example). If the upstream and downstream velocities 
are equal, their value and the value of some other velocity 
(the maximum or minimum welocity, for example) can be 
selected for q a and q b ; or, if desired, q a can be equal to q b , in 
which case if 

q a —q-\-e where e -»0 
2>=2 


it can be shown from equations (14a) and (14b) that 


and 



(C4a) 


(C4b) 



From equation (13) and the conditions given by equation (C3) 



This latter case, in which q a =qb=q, corresponds to the 
method used by Chaplygin (ref. 6 ) and Kdrmdn-Tsien (ref. 8 ) 

in which the correct relation between p and — is replaced 

p 

by a straight line (eq. (Cla)) that is tangent to the correct 
relation at one point (where q a =qt)- 


APPENDIX D 

EQUATIONS OF CONTINUITY AND ERROTATIONAL FLUID MOTION IN TERMS OF TRANSFORMED <p*, ^-COORDINATES 


Consider the two-dimensional irrotational motion of a 
fluid particle in the physical xy-plane. The fluid particle is 
defined by adjacent streamlines (constant 4'*) and velocity- 
potential lines (constant <p*) spaced 5n and 5s apart as indi- 
cated in figure 23. The velocity q* is parallel to the stream- 
lines and normal to the velocity-potential lines. 

Continuity. — From continuity considerations of the fluid 
particle in figure 23 


b log, p* d<p* . d log, g* d<p* . d0 df* 
dip* ds d^* ds dn 


or, from equations (15) and (16) 

1 / d log, p* d log, q* \ . d0 

^ V V d<p* J W 


=0 


But, from equation (13) 


(Dl) 


A (p V5n )= 0 

or 

d log, p* d loga q* . 1 b(Sn)_ 
ds ds ds 


1 d log, p* _ — q* s dlog, q* 
P* dp* Vl + 2** 

so that equation (Dl) becomes 

1 d log, g.* be 

Vl + g 55 V- W 


which combined with equation (B2a) becomes 


(D2) 
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U= 3—. 

l+VT+s 55 

then 

(18) 

d>3) 

so that equation (D2) becomes 


b log, u , be n 
b P * + bip*~° 

(17) 

Equation (17) is the continuity equation expressed in terms 
of <p*, ^-coordinates and log, u. 

Irrotational fluid motion. — For irrotational motion of 
the fluid particle in figure 23 




or 


a, log, g* . 1 b(0s) ^ Q ~ 
bn bn 

which, combined with equation (B2b) becomes 

b log, q* dip* b6 d<p* q 
bp* dn b<p* ds ~ 

or, from equations (13), (15), and (16) 

1 . b log, q* b0 _ n 

■ VT+?^ W 


Finally, from equations (D3) and (D4) 
b log, u b9 




btc* 


0 


(D4) 

( 20 ) 


Equation (20) is the equation for irrotational fluid motion 
expressed in terms of <p*, ^-coordinates and log, u. 


APPENDIX E 


INTEGRAL EQUATION FOR 0($„*,>) 


If the distribution of the angle 0(4?, 1 $ r ) in the transformed 
d'SF-plane is harmonic, that is, satisfies equation (45) within 

and on the channel walls (*St equals 0 and then from 

Green’s theorem and the theorem of mean value it can he 
shown that the value of 6 at a point (4 ) 0 ,'4 r <) ) within (or on) 
the channel walls is given by (ref. 9, p. 204, for example) 



where the two integrals on the right side of equation (El) 
represent the line integral around the channel walls in the 
counterclockwise direction with the signs adjusted so that 


^ represents the inner normal to the path of integration. 


The function £?(<f>,'4 r ) in equation (El) is of the form 
(ref. 9, p. 204) 

<?($,*) =log, -+ cu ($,¥) (E2) 

T 


where r is the distance from any point (4?,’3i r ) to the point 
('f> 0 ,4'' 0 ) and where co($,4') is an arbitrary function that is 
harmonic within and on the channel walls. (Thus from 
equation (E2), .6r(4>,4 r ) is harmonic within and on the channel 
walls except at the point (4> 0 , 1 4 r 0 ) where a logari thmi c singu- 
larity exists.) Because the harmonic function <o(4>,4 r ) is 
arbitrary, the function (?(4>,'$ r ) can be selected so that along 


the channel-wall boundaries 


(* equals 0 and 0 || i 


is a con- 


stant c given by the following equation (obtained from notes 
presented by Tamarkin and Feller in the 1941 Summer 
Session for Advanced Instruction and Research in Mechanics 
at Brown Univ.) : 



(E3) 


where l is the length of the path along which the line integral 
is taken. For the path under consideration l is infinite and 

therefore (7(4?, 4r) can be selected so that ^ is zero along the 

channel walls. A function with this property is a Green’s 
function of the second kind. Equation (El) becomes 


0(4?o,*«)=o- 


L 2 




or, combined with equation (43a) 


/: r(« ^xi * <«> 

La 


Along the channel walls — is known from the prescribed 

velocity distribution so that, after the proper Green’s func- 
tion O has been determined (appendix F), equation (46) de- 
termines the value of 9 at any point ($ £ ,,4 r ,). The valuo of 
0(4> o ,^o) given by equation (46) can be adjusted by an arbi- 
trary constant of integration to give a specified value of 0 at 
one point in the flow field. 
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APPENDIX F 


GREEN’S FUNCTION OF SECOND KIND 


From appendix E Green’s function of the second kind Q 
satisfies the condition 

^=0 
d'f r U 

along the channel walls, which are straight and parallel 
boundaries Ak equals 0 and ^ extending to ± °°‘in the <f>- 
direction, and satisfies the equation 

d 2 (? , d 2 g 
&<k 2 + 

everywhere in the channel except at the point (# 0 ,"k 0 ) 
where Q has a logarithmic pole. For these conditions the 
Green’s function Q can be obtained by analogy from the 
velocity potential for incompressible flow into a point sink at 
(ko/k,) between straight parallel boundaries at k equal to 0 

and ~ The logarithmic pole for Q at (k^k,) corresponds 
z 

to the point sink, and the condition ^=0 at the boundaries 

corresponds to zero velocity, that is, no flow normal to the 
boundaries. 

The velocity potential for fluid flow with the boundary 
conditions just described is obtained from two infinite 
series of point sinks with the sinks of each series spaced * 
distance apart in the ^-direction, and the two series arranged 
by the method of images in such a manner that no flow crosses 

dG 

the boundaries, that is, ^-^=0. This arrangement of point 
sinks is shown in figure 24. 

The complex function Wi for the first infinite series of 
point sinks is given by (ref. 10, p. 112, for example) 


where 


Wi=— log, sinh ( z—z „ ) (Fla) 

z=k-Hk (Fib) 

The complex function Wj for the second infinite series of 
point sinks ^mirror image of the first series in order to prevent 

flow across the boundaries k equals 0 and ^ is given by 

Wi=— log, sinh (z—z 0 ) (F2a) 

i=k— ik (F2b) 


where 


The complex function w for the combined flow becomes from 
equations (Fla) to (F2b) 

w=vh+ith=— log, sinh [(<!>— k<,)-H(k—k,,)]— 

log, sinh [($—$,) -H(k+k„)] (F3) 


Point sinlts given by 



Fiqctee 24.— Two infinite series of point sinks required in the development of Green’s function 
of tho second kind O. 

The Green’s function of the second kind 6 corresponds to 
the velocity potential for the incompressible flow and is there- 
fore given by the real part of equation (F3) 

Q=— — log, [cosh 2 (k— $„)— cos 2 (k— k 0 )][cosh 2 ($—$,)— 

cos 2 (k+k„)] (F4) 

But along the channel walls k is equal to 0 or ~ so that 

z 

cos 2 (k+k 0 )=cos 2 (k— k 0 ) 
and equation (F4) becomes 

O 0or t=— log, [cosh 2 (k— k„)— cos 2 Ok— 4',)] (47) 

Equation (47) gives the Green’s function of the second kind 
along the channel walls (straight parallel lines of constant k 

equal to 0 and ^ and extending to ± m in the k-direction) . 
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Several techniques, depending on the magnitude of the 
upper limit |($— <£,)|, were used to evaluate the integrals 
a and P given by equations (50d) and (50e). Each integral' 
is treated separately in this appendix, and the values of 
($—$„) f or the upper limit [($—$„)] are considered positive. 
For negative values of (<£—$„) the magnitudes of I (that 
is, of a or p) are equal for corresponding values of | ($—$„) | 
but opposite in sign. As a result the values of Al have the 
same sign. 

. INTEGRAL a 

Small and medium values of (4> — $ 0 ). — For small and 
medium values of the upper limit of integration ($— $„) 
in equation (50d), that is, for 0 <($>—$„)< 60*724, the 
integral a is evaluated by Simpson’s one-third rule using 
increments of ($—$„) equal to x/48. 

large values of -(4>— <&„).. — For large values of ($—$,), 
that is, for ($— $o)>60+24, the integrand in equa- 
tion (50d) becomes 


APPENDIX G 

EVALUATION OF a AND 0 
(*-*,) 


log, cosh log, 2 

so that equation (50d) becomes 

/’Mr/34 ‘ 

a« I log, cosh ($—$,) </($—$,)+ 
f ( 0> [($—$,)— log, 2] </($—$„) 

* i/6 Or/34 


(Gl) 


*25.809782+ 


»0.411230-0.693147(<i»-$ ( ,)+£ ($-$,)* (G2) 


-0.693 147(4> — $,) — 25.398552 


0 


1 

2 


Equation (G2) gives values of a for values of ($—$„) 
equal to or greater than 60+24. Values of the integral a are 
tabulated in table VJJL for a range of | ($ — $,) | between 0 
and 100+24 in increments of +24. For negative values 
of ($—$,) the sign of a is negative. 

INTEGRAL 0 

Small values- of ($ — 4*<>). — For ($ — $,) equal to zero the 
integrand of equation (50e) becomes infinite so that Simpson’s 
one-third rule cannot he used to evaluate & in this region of 
($—$,), as was done for a. However, equation (50e) 
integrates by parts to give 


x: 


log, sinh ($—$„) d(4>— $,)=($— 4>,) log, sinh (4>— <!>,)— 


X 




($—$,) ctnh ($—$„) d(4>— $„) (G3) 


where the integrand ($—$,) ctnh ($— 4> 0 ) on the right side 
of equation (G3) can be expanded in the following series form : 

(*-*■) ctah <#-*.) 


2 6 H 5 (4>- $,) 9 2 8 H 7 ($-$,) 8 2 10 B 0 (4>— $,) 10 
6! 8! + 10! 

2 12 B U (4> — $ # ) 12 
12 ! 


... (G4) 


where B 1} B z , and so forth, are Bernoulli’s numbers (ref. 11, 
p. 90, for example). From equations (G3) and (G4) 


j3 = ($— $,) log, sinh ($-$„)—($—$„)- 


($-$ 0 ) 3 . (4-4>,) 5 


-. . . (G5) 


9 1 225 

2(4> — $ 0 ) 7 (3>-3>,) 9 2(^-$ 0 ) n 1382(4— I-,) 13 

6615 '42,525 l,029,105 + 8,300,667,375' 

Equation (G5) was used to obtain p as a function of ($—$„) 
for 0 <($—*„) <8+24. 

Medium values of ($— $„). — For medium values of the 
upper limit of integration (4>— <£„) in equation (50e), that is, 
for 8+24<($— $„)< 60+24, the integral p is evaluated by 
Simpson’s one-third rule as was done for a. 

Large values of ($—$„). — For large values of ($—$,), that 
is, for ($—$,) >60+24, the integrand in equation (50e) 
becomes 

log, sirih ($+$„) « ($—$,) —log, 2 (G6) 

so that equation (50e) becomes 

r /24 

log, sinh (<£—$„) +4>— +,)+ 


» 

Ja 


(*- *„) 


60x/24 


[($—+)— log, 2] d ($— 4>„) 


«24.576082+p-^- ! -0.693147(4>-4\,)-26.398562j 

«-0.822470-0.693147(4>-$,)+i ($-$,)* (G7) 

Equation (G7) gives values of p for values of (4>— 4>,) equal to 
or greater than 60+24. Values of the integral p are tabulated 
in table VII for a range of )($— 4>,)| between 0 and 100+24 
in increments of +24. For negative values of (4>— 4?,), the 
sign of p changes. 


APPENDIX H 

CHANNEL TURNING ANGLE 


If the prescribed velocity distribution along one channel 
wall differs from the distribution along the other wall, then 
in general the channel deflects an amount A 6, which is the 
difference in flow direction far downs treajn and far upstream 


of the region in which the prescribed velocity distribution 
varies. Thus, 


A0=6i — 6-u 


(HI) 
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CH, 


For large values of | (4>— $ 0 ) | such as occur far upstream 
and far downstream of the region in. which the prescribed 
velocity varies along the channel walls 


(H2) 


cosh 1 2 3 4 * 6 ($—$„)>> cos 2 

SO that from equation (47) 

/ G 0 =Or= — 2[|($— $ 5 )|— log, 2] 

2 

Far upstream so that 

| ($-$.) | = (*-*„) 

and because V is harmonic 


so that equation (H2) substituted into equation (46) gives 


Likewise, far downstream so that 

| ($-$„)[ = -($-$„) 

and equation (H2) substituted into equation (46) gives 

2 

From equations (Hi), (H3), and (H4) 

2 

Equation (H5) determines the ch ann el turning angle A 9. 
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TABLE I— DISTRIBUTION OF VELOCITY Q AND FLOW DIRECTION B IN TRANSFORMED ^t-PLANE FOR EXAMPLE HI 

(ELBOW WITH INCOMPRESSIBLE FLOW) 

[Prescribed variation in Q with aro length t along channel walls plotted in flg. 2; <J»— 0.6, 0 j**1 J>, A«-89.33°] 
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TABLE II— DISTRIBUTION OF PHYSICAL COORDINATES a AND y IN TRANSFORMED w^PLANE FOR EXAMPLE in (ELBOW 

WITH INCOMPRESSIBLE FLOW) 

[Prescribed variation in <? with aro length t along channel walls plotted In fig. 2; 0.— 0.5, <?j=>1.0, Afl-.S9.3S 0 ] 
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0 
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0.250 
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0.760 


1.000 

w \ 

X 

V 

X 

V 

X 

V 

D 

V 

n 

D 

B 

D 

D 

D 

D 

D 

B 

KM 

-2.000 
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-a 078 

-a 748 

-a 978 

-a 498 

-2978 

-a 248 

-3 978 

a 002 

-3978 

0.252 

-3 978 
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a 752 
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-3.727 
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-2 723 
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-3 728 
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.252 
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-2478 
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-a 478 
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.252 

-3478 
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1.002 
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-3 227 

—.098 
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BH 
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0 
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.000 
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.000 

.000 
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.125 
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.254 
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.243 

-.476 

.239 

-.239 

.240 

-.002 
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.252 
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.672 
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-.254 
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-.041 
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843 
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.607 
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-.031 
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-.044 
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-.316 
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.055 
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—.849 
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2422 

-.658 

2.504 

-.398 
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-.228 
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-.044 
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a 413 
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2488 
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a 768 
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.030 
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2.535 
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-.594 

3930 
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3056 

-.281 
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1.875 
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2741 

-.900 


-.838 

2.957 
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-.418 
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2964 
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3501 
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2.250 
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2172 
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3.465 

-1.002 
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am 
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3476 
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TABLE HI — DISTRIBUTION OF VELOCITY q AND FLOW DIRECTION 8 IN TRANSFORMED ptyAPLANE 
FOR EXAMPLE IV (ELBOW WITH LINEARIZED COMPRESSIBLE FLOW) 

(Prescribed variation in Q with arc length s along channel ’■rails plotted In fig. 2; <J.— 0.6, Qd- '1-0, ji—0.80176, 0.73782, AS= 104.07°) 
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TABLE IV— DISTRIBUTION OP PHYSICAL COORDINATES * AND y IN TRANSFORMED p*^*-PLANE FOR EXAMPLE IV 

(ELBOW WITH LINERARIZED COMPRESSIBLE FLOW) 

(Prescribed variation in <? with arc lengthsalong channel wallsplottedinflg. S; Q.-0.5, Oj— 1.0, gj=0^0176, 0.73782, A0-1O4.O7 0 ] 
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TABLE V— DISTRIBUTION OF VELOCITY q AND FLOW DIRECTION 6 IN TRANSFORMED p^PLANE FOR EXAMPLE V (ELBOW 

WITH COMPRESSIBLE FLOW (y =1.4)) ' . 


[Prescribed variation In Q with arc length < along channel walls plotted in fig . 2 ; <?.— 0 . 6 , Q .- LO , J .- 0 . 79927 , Af - 0 . 71054 , Ai - 105 . 31 0 ] 


\ * 

0 

H 

H 

a 



1.0 

<P \ 

ff 

0 

Q 

e 

? 

e 

$ 

9 

1 

6 

ff 

6 


9 

— 12/6 

CL 3990 

0 

a3996 

0 

a 3990 

0 

a 3990 

0 

03996 

0 

0 3990 

0 

HI 

0 

— 11/6 

.3990 

.00 

.3997 

.00 

.3997 

.00 

.3997 

.00 

.3997 

.00 

.3997 

.00 


.00 

- 10/6 

.3996 

.01 

.3997 

.01 

.3997 

.00 

.3997 

.00 

.3997 

.00 

.3997 

iim 

.3990 

-.01 

- 9/6 

.3990 

.01 

.3997 

.01 

.3997 

.01 

.3997 

.00 

.3997 

-.01 

.3997 



1 iffiOr/W 

-are 

.3990 

.02 

.3997 

.02 

.3998 

.01 

.3998 

.00 

.3998 

-.01 

.3997 


.3990 


- 7/0 

.3996 

.03 

.3998 

.03 

.3999 

.02 

.3999 

.00 

.3998 

-.02 

.3998 

-.03 

.3996 


-m 

.3990 

.05 

.3998 

.01 

.4000 

.03 

.4000 

.00 

.4000 

-.03 

.3998 


.3990 


- 6/0 

.3990 

.09 

.4000 

.08 

.4002 

.04 

.4003 

.00 

1 .4002 

-.04 

.4000 

^fcrra 

.3996 

-.08 

— 1/6 

.3996 

.15 

.4002 

.13 

.4006 

.07 

.4008 

.00 

.4006 

-.07 

.4002 

-.12 

.3998 

-.14 

- 3/6 

.3996 

.25 

.4007 

.22 

.4014 

.12 

.4015 

—.01 

.4012 

-.12 

.4005 

-.20 

.3996 

-.23 

- 2/0 

.3990 

.43 

.4015 

.30 

.4020 

.18 

.4028 

-.03 

.4022 

—.22 

.4011 

-.34 

.3990 

-.38 

- 1/8 

.3996 

.77 

.4030 

.62 

.4047 

.30 

.4019 

-.06 

.4038 

-.36 

.4019 

-.65 



0 

.3996 

1.46 

.4062 

1.06 

.4080 

.39 

.4032 

-.19 

.4062 

-.64 

.4031 

-.92 

.3998 

- 1.02 

1/8 

.4060 

a 13 

.4137 

1.01 

.4152 

.40 

.4134 

-.45 

.4097 


.4049 

- 1.48 


- 1.01 

2/6 

.4263 

4.28 

.4276 

. 1.93 

.4255 

.24 

.4207 

—.98 

.4144 

- 1.81 

.4072 

— 2.30 

.3990 

- 2.40 

3/6 

.-4519 

4.31 

.4461 

1. 54 

.4389 

-.48 

.4300 

- 1.92 

.4202 

- 2.89 

.4099 

-046 

.3990 


4/6 

.4830 

a 20 

.4687 

.40 

.4517 

- 1.74 

.4408 

—0 30 

.4268 

- 4.35 

.4131 

- 4.97 

l ! MVcTL ! Ti^H 

- 5.17 

6/6 

.6102 

1.21 

.4928 

- 1.53 

.4719 

-a 04 

.4524 

-0 21 

.4340 

-0 29 

.4104 

-0 92 


- 7.13 

0/6 

.6499 

- 1.01 

.5175 

- 4.11 

.4895 


.4013 

- 7.61 

.4413 

- 8.60 

.4198 

- 9.28 


- 0.49 

7/6 

.5823 

-a 05 

.6119 

- 7.27 

.6069 


.4762 

- 10.50 

.4485 

- 11.49 

.4231 



- 12.28 

8/6 

.0144 

-a 99 

.5652 

- HL92 

.5236 

— 12.55 

.4875 

-1082 

.4654 

- 14.74 

.4263 

-1029 


-1047 

9/6 

.0141 

— 13. 32 

.6871 

- 14.98 

.6393 

- 16.40 

.4981 

- 17.51 

.4618 

— 18.37 

.4292 

— 18.80 

.3990 


10/6 

. 6717 

- 17.98 

.6074 

- 19.38 

.6538 

- 2a 61 

.6078 

- 21.61 

.4677 

- 22.34 

.4319 

- 22.79 

.3998 

- 22.91 

11/0 

.6970 

- 22.89 

.6258 

- 24.06 

.5669 

— 25. 12 

.5166 

-20 98 

.4730 

-20 02 

.4343 


.3990 

- 27.16 

12(6 

.7197 

- 28.02 

.6424 

- 28.98 

.6788 

- 29.87 

.6246 

- 3a 62 

.4777 

- 3L18 

.4364 

- 31.62 

.3990 

- 31.03 

13/6 

.7399 

- 3a 31 

.6569 

— 34. 09 


- 34.84 

.5314 

-3048 

.4818 

-3096 

.4383 

-30 28 

.3990 


14/6 

.7673 

- 38.74 

.6690 

- 39.35 

.6978 

- 39.98 

.5376 

— 4a 52 

.4855 

- 4a 94 

.4400 

- 41.20 

.3990 

- 41.30 

16/0 

.7719 

- 44.28 

.6803 

- 44.74 

.6056 

— 4a 25 

. 5428 

-40 72 

.4888 

-40 09 

.4415 

- 46.33 


-4541 

ie/6 

.7836 

— 49.&1 

.6891 


.6123 

- 5 a 62 

.6170 

- 5L03 

.4918 

- 51.37 

.4429 

- 61. 61 


- 61.70 

17/0 

.7922 

- 65.44 

.6963 

- 55.70 

.6182 

- 5a 04 

.5522 

-6043 

.4951 

-5080 

.4446 

- 57.07 


- 67.18 

18/0 

.7975 

— 00.96 

.7019 

- 61.13 

.6238 

— 61.44 

.5573 

— 61,85 

.4990 

- 62.30 

.4468 

- 62.69 

.3998 

- 02.80 

10/0 

.7983 

- 60.35 

.7066 

— 66.45 

.6297 

-00 76 

.5635 

- 67.25 

.5047 

- 67.87 

.4607 

- 68.67 


- 08.97 

20/6 

.7993 

- 71.43 

.7110 

- 71.53 

.6307 

- 7L8S 

.5720 

- 72.51 

.5138 

- 7a 42 

.4695 

- 74.59 



21/6 

.7993 

- 70. 19 

.7163 

- 70.31 

.6456 

- 7a 73 

. 6835 

- 77.60 

.5274 

- 78.69 

.4753 

- 80.35 

.4253 


22/6 

.7993 

- 8a 67 

.7224 

it. ■■ 

.6562 

- 81.27 

.6979 

- 82.14 

.5453 

-8048 

.4009 

-8046 

.4619 

- 88.20 

23/0 

.7993 

— 81.69 

.7292 

— 84.83 

.6684 

—8033 

.6145 

-80 27 

.5662 

- 87.72 

.5225 

- 89.70 

.4830 


24/0 

.7993 

- 83.32 

.7366 

- 88.47 

.6816 

- 88.97 

.6327 

- 89.92 

.5892 

- 91.38 

.5603 

- 03.41 

.5102 

-9511 | 

26/6 

.7993 

- 91.56 

.7442 

- 91.09 

.6952 

- 92.17 

.6616 

-90 06 . 

.6129 

- 94.41 

' .'5789 

-90 29 

. 6199 

- 98.75 

26/0 

.7993 

- 91.33 

.7516 

- 91.46 

.7087 

- 94.90 

.0705 

- 95.72 

.6366 

-90 96 

.6073 

- 98.05 



27/8 

.7993 

- 06.70 

.7687 

— 86 l 81 

.7219 

- 97.21 


- 97.94 

.6597 

- 99.04 

. 0*47 

- 100.62 

.0144 

- 102.30 

28/6 

.7993 

- 98.68 

.7853 

- 98.78 

.7342 

- 99.14 

.7062 

- 99.77 

.6316 

-loan 

.6607 

- 101.90 

.0141 

- 103.57 

29/0 

.7993 

- 100.31 

.7714 

-10940 

.7456 

-loan 

.7223 

- 101.25 

.7019 

- 102.00 

.6849 

- 103 . J3 

.0717 

- 101.46 

30/6 

.7993 

- 101.01 

.7709 

- 101. 72 

.7559 

- 101.97 

.7370 

- 102.43 

-.7205 

-10010 

.7070 

- 103.99 

.0970 


31/6 

.7993 

- 102 . 6S 

.7817 

- 102.75 

.7651 

- 102.97 

.7502 

- 103.34 

.7373 

- 103.89 

.7270 

- 101.01 

.7197 

- 105.40 

32/0 

.7993 

- 103.48 

.7858 

- ICC . 51 

. .7731 

— 100 72 

.7617 

- 104.03 

.7520 

- 104.47 

.7446 


.7399 

-105 76 

33/6 

.7993 

- 101.09 

.7893 

- 101.14 

.7800 

— 104.28 

.7716 

- 10L63 

.7048 

- 101.88 

.7699 

- 105.33 

.7673 

-10587 

31/6 

.7993 

— 104.63 

.7922 

- 104.56 

.7856 

- 10L67 

.7799 

- 104.87 

.7754 

-10014 

.7726 

-10049 

.7710 

- 105.89 

36/0 

.7993 

— 101.83 

.7915 

- 101.86 

.7901 

- 10L95 

.7865 

-100 09 

.7839 

-10029 

.7828 

- 105.61 

.7830 

-10581 

36/8 

.7993 

- 105.04 

.7962 


.7935 

-10011 

.7914 

- 106.21 

.7903 

-10035 

.7904 

-100 52 

.7922 

-100 71 

37/6 

.7993 

- 105.16 

.7975 

- 106.17 

,.7959 


.7948 

-100 28 

.7945 

-10036 

.7953 

-10046 

.7076 

- 105.50 

38/6 

.7993 

- 105.22 

.7983 

- 105.21 

.7974 

-10026 

.7969 

-10030 

.7970 

-10034 

.7977 

- 105.39 

.7993 

-100 41 

39/6 

.7993 

- 105.20 

.7887 

- 105.27 

.7983 

-100 28 

.7981 

-100 31 

.7982 

-10033 

.7986 

-10035 

.7993 

-10035 

40/6 

.7993 

- 105.29 

.7990 

- 105.29 

.7988 


.7987 

-10031 

.7987 

-10032 

.7990 

- 105.33 

.7993 

- 105.33 

41/0 

.7993 

- 105.30 

.7991 

- 105.30 

.7990 

E iL f'l 

.7990 

-100 31 

.7990 

-100 31 

.7991 

-100 32 

.7993 


42/6 

.7993 

- 105.31 

.7992 

- 105.31 

' .7992 

-100 31 

.7991 

- 105.31 

.7992 

-10031 

.7992 


.7093 


43/0 

.7993 

- 105.31 

.7993 

- 105.31 

.7992 

-100 31 

• .7992 

-10031 

.7992 

-10031 

.7993 


.7993 

- 105.31 

44/6 

.7993 

- 105.31 

.7993 

- 105.31 

.7993 

-10031 

.7993 

-10031 

.7993 

-10031 

.7993 

-10031 

.7993 

-10031 
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TABLE VI— DISTRIBUTION OF PHYSICAL COORDINATES x AND y IN TRANSFORMED W--PLANE FOR EXAMPLE V (ELBOW 

WITH COMPRESSIBLE FLOW ( 7 =1.4)) 

[Prescribed variation in 0 with arc length j along channel walls plotted in flg. 2; 0.— 0.6, Os— 1.0, gso 0.79927, A^-0.71054, AS=106^1°] 
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TABLE VII— TABULATED VALUES OF THE INTEGRALS a AND 0 FOR A RANGE OF |($— $„)| 

[Computational methods given In appendix G] • 


1(4— $.)| 

aW 

0 

0 

l(x/24) 

.000373 

■ 2(x/24) 

.002970 

3(x/«) 

.009942 

4(x/24) 

.023301 

5(x/24) 

.044875 

6(x/24) 

.070277 

7(x/24) 

.118897 

8(x/24) 

.173909 

9(x/24) 

.242283 

10(x/24) 

.324812 

ll(x/24) 

.422136 

12(x/24) 

.534763 

13(x/24) 

.'66309S 

14(x/24) 

.807460 

16(x/34) 

.968096 

lfl(x/24) 

L 146201 

17(x/24) 

L 338926 

lS(x/24) 

L 549389 

19(x/24) 

L 776679 

20(x/24) 

2.020867 

21(x/24) 

2.282008 

22(x/34) 

2.560143 

23(x/24) 

2.865304 

. 24(x/24) 

1167616 

25(x/34) 

3. *90799 

28(x/24) 

3.843168 

27(x/24) 

4.206633 

28(x/24) 

4.587203 

29(x/24) 

4.984835 

30(x/24) 

5.399685 

31(x/24) 

5.831607 

32 (x/24) 

6.280652 

33(x/24) 

6.746825 

34(x/24) 

7.230126 

35(x/24) 

7.730557 

38(x/24) 

8.248119 

37(x/24) 

8.782813 

33(x/24) 

9.334640 

39(x/24) 

9.903600 

40(x/24) 

la 489693 

41(x/24) 

11.092920 

42(x/24) 

1L 713281 

43(x/24) 

12.350777 

44(x/34) 

13.005406 

4S(x/24) 

13. 677170 

46{x/24) 

14.366063 

47 (x/24) 

15.072101 

48(x/24) 

16.795268 

49 (x/24) 

16.635571 




Al»Arz 


Aia 

(Afr*«724) 


a 000373 


.002697 


.006972 


.013369 


.021574 


.031402 


.042620 


.055012 


.068374 


.082529 


.097324 


.128336 


.144362 


.160636 


.177105 


.193726 


.210463 


.227290 


.244188 


.261141 


.278136 


.295161 


.312212 


.329283 


.346369 


.363466 


.380570 


-.397682 


.414800 


.431922 


.449045 


.466173 


.483301 


.600431 


.517562 


.634694 


.651827 


.568960 


.586093 


.603227 


.620361 


.637496 


.654629 


.671764 


.706033 


.723167 


.740303 


.757438 


A ter 

(AS-2r/24) 


0.002970 


.020331 


.052976 


.097632 


.150903 


.209951 


.272697 


.337741 


.404188 


.471478 


.539276 


.607373 


.675652 


.744035 


.812482 


.880967 


L 017993 


L 086521 


L 155053 


L223588 


L 292125 


L 380682 


1.429200 


L 497739 


£U) 

0 

-.396937 

— . 611660 

-.756406 

-.854451 

-.916486 

-.948630 

-.954346 

-.936438 

-.896542 

-.836896 

-.755399 

— . 655715 

-.637338 

-.400640 

-.245901 

-.073335 

.116897 

.324671 

.549892 

.792488 

1. 052403 

L 329594 • 

L 624029 

1.935683 

X 264536 

2.610573 

2.973783 

3.354158 

3.751689 

4.166373 

4.598205 

5.047181 

5.513301 

6.996561 

6.490961 

7.014499 

7. 549176 

8.100987 

8.669936 

9.256021 

9.859241 

10.479597 

11. 117089 

1L 771715 

12.443477 

13.132373 

13.838405 

14.561571 

15.301873 



AI-A0 


Ai/1 

(A*- */24) 

Aip 

(A4-- 2t/24) 

-a 396937 

-a 611660 

-.214723 

144746 

-.242791 

-.098045 

-.062036 

-.094079 

— . 032044 

-.006816 

.012092 

.017908 

.039896 

.100542 

.060640 

.080497 

.180181 

.099684 

.118377 

.255075 

.136698 

. 154739 

.327305 

.172566 

.190232 

.398006. 

.207774 

.'225221 

.467817 

.242596 

.259915 

.637106 

.277191 

.294435 

.606089 

.311654 

.328853 

.674890 

.346037 

.•363210 

.743585 

.380375 

.397631 

.812216 

.414684 

.431832 

.880808 

.448976 

.466120 

.949380 

.483260 

.600100 

L 017938 

.617538 

.534676 

1.086483 

.551812 

.568949 

1.156034 

: 686085 

.003220 

L 223576 

.620366 

.637492 

1*292118 

654626 

.671762 

L 360658 

.688890 

.706032 

L 429198 

.723166 

.740302 

h 497738 

.757436 


to For negative values of (<£— <£•) the signs of a and p change, but the signs of Aa and remain unchanged. 
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TABLE VTI— TABULATED VALUES OF THE INTEGRALS a AND p FOR A RANGE OF |(«-$„)|— Concluded. 

[Computation methods given in appendix G.] 
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192 REPORT 1115 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 

TABLE VIII— COMPARISON OF ELBOW DESIGNS OBTAINED FROM SOLUTIONS BY RELAXATION METHODS AND BY 

GREEN’S FUNCTION 

[Linear lied compressible flow; prescribed velocity distribution given In figs. 2 and 22.] 
























